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We construct a theory for Bose-Einstein condensation of light in nano-fabricated semiconductor 
microcavities. We model the semiconductor by one conduction and one valence band which consist 
of electrons and holes that interact via a Coulomb interaction. Moreover, we incorporate screening 
effects by using a contact interaction with the scattering length for a Yukawa potential and describe 
in this manner the crossover from exciton gas to electron-hole plasma as we increase the excitation 
level of the semiconductor. We then show that the dynamics of the light in the microcavities 
is damped due to the coupling to the semiconductor. Furthermore, we demonstrate that on the 
electron-hole plasma side of the crossover, which is relevant for the Bose-Einstein condensation of 
light, this damping can be described by a single dimensionless damping parameter that depends 
on the external pumping. Hereafter, we propose to probe the superfluidity of light in these nano- 
fabricated semiconductor microcavities by making use of the differences in the response in the normal 
or superfluid phase to a sudden rotation of the trap. In particular, we determine frequencies and 
damping of the scissors modes that are excited in this manner. Moreover, we show that a distinct 
signature of the dynamical Casimir effect can be observed in the density-density correlations of the 
excited light fluid. 


I. INTRODUCTION 

The first observation of Bose-Einstein condensation in 
a dilute atomic vapor opened up several different possibil¬ 
ities to explore many-body phenomena in a completely 
new regime. As before this observation it was hard to 
experimentally access the macroscopic quantum regime, 
the low temperatures of these systems created a play¬ 
ground for the investigation of several interesting quan¬ 
tum effects. One prominent example is the observation 
of superfluidity via the existence of quantized vortices 
m- These developments were even more encouraged 
by the large experimental control that is available in the 
cold atomic gases, leading also to the possibility to inves¬ 
tigate dynamical behavior. 

More recently a class of Bose-Einstein condensates 
have been observed that are fundamentally different 
from atomic condensates and also allow for different 
experimental probes. These Bose-Einstein condensates 
of bosonic quasiparticles such as magnons [1], exciton- 
polaritons and photons [7] are dissipative systems. 
This leads to the possibility to investigate phenomena 
that are not not yet observed in atomic condensates, such 
as large number fluctuations in a Bose-Einstein conden¬ 
sate of photons mm- Moreover, it is also interesting to 
investigate whether certain equilibrium phenomena are 
still present in these dissipative systems. 

Although the Bose-Einstein condensate of light al¬ 
lows for dynamical measurements and new experimental 
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probes, there are still some disadvantages to the current 
experimental approach [3 HU]. First, it is inconvenient 
to dynamically change the trap geometry of the photons. 
Second, every possible change highly affects the interac¬ 
tions of the photons with the dye and therefore strongly 
influences the thermalization of the photons. Lastly, it is 
very hard to predict theoretically, and to a large extent 
it is still unknown experimentally, how strongly the pho¬ 
tons interact and what the origin of the effective photon- 
photon interaction is. Therefore, this is a large disad¬ 
vantage for studying various many-body phenomena of 
light. 

To overcome these problems we propose to follow a 
different experimental approach, namely nano-fabricated 
semiconductor microcavities. In this system we use a 
semiconductor with a periodic array of holes filled with 
air to create a photonic crystal. The photons thermalize 
due to the interaction with the electrons and holes in the 
semiconductor. Moreover, by simply increasing the size 
of the holes when moving further away from the center of 
the semiconductor, the photons feel an effective harmonic 
trapping potential m- Finally, also in this case there is 
external pumping to compensate for photon losses out of 
the microcavity. Although this suggest that light in such 
nano-fabricated semiconductor microcavities is similar to 
photons in a dye-filled optical microcavity, the former 
system has several advantages. Most importantly, the 
interaction between the light and the semiconductor is 
well understood, see e.g. Ref. [TU]. Therefore, there are 
good prospects of controlling the effective photon-photon 
interaction. Moreover, it is possible to change the cavity 
geometry and it is also achievable to influence the photon 
trap dynamically on a sub-picosecond timescale [T5] . 
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A prime example of a many-body phenomenon of 
light where the advantages of the semiconductor micro¬ 
cavity can be utilized is superfluidity. In Ref. [T3] it is 
shown that in dissipative condensates a superfluid den¬ 
sity can exist and therefore it is expected that also a 
Bose-Einstein condensate of photons can exhibit super¬ 
fluidity. Although some work has been carried out that 
mainly focus on the frictionless flow of light through an 
obstacle, see e.g. Ref. [TMT^ . up to now the superfluid¬ 
ity of a Bose-Einstein condensate of photons has not yet 
been explored, since it is unclear how the superfluidity 
of this system can be probed in the current experiments. 
Namely, recall that the standard experiment for probing 
the superfluid behavior, i.e., observing the existence of 
quantized vortices after rotation of the condensate, ap¬ 
pears not to be feasible in a dye-filled optical microcavity. 

From studies in atomic Bose-Einstein condensates we 
know that there is another method to obtain direct evi¬ 
dence for superfluidity [IIIID]. By applying a small sud¬ 
den rotation of the trapping potential the so-called scis¬ 
sors modes can be excited. The modes are quite different 
from the excitations in the of the atoms in the normal 
state. Therefore, the time evolution of the angle between 
the axial direction of the condensate and the new trap di¬ 
rection is distinct in the normal and superfluid phase. As 
mentioned before, in current experiments this is difficult 
to observe as a rotation of the cavity highly affects the 
interactions of the photons with the dye and therefore 
destroys the thermalization of the photons. In contrast, 
in a semiconductor microcavity the trap geometry can 
be changed on a sub-picosecond timescale [T3], which is 
much shorter than the inverse trapping frequency which 
is typically on the order of tens of gigahertz. Therefore, 
this fast change of the cavity geometry is non-adiabatic, 
which allows for the excitation of the scissors modes and 
probing the superfluidity of a Bose-Einstein condensate 
of light. 

Nonetheless, there are still major differences between 
the dynamics of the scissors modes in Bose-Einstein con¬ 
densates of atoms and photons. The former systems are 
very clean and therefore damping is typically not an is¬ 
sue, especially because most experiments are carried out 
in the collisionless regime. As a result, only for fine- 
tuned configurations of the experimental set-up, such as 
in Ref. [H], the scissors modes are damped by Beliaev 
processes. In a Bose-Einstein condensate of light on the 
other hand, the photons are, as envisaged here, coupled 
to electrons and holes in the solid-state cavity, that is 
pumped with an external laser beam. In the context 
of exciton-polariton systems it is shown that this pump¬ 
ing of the external bath affects the superfiuid properties 
[221 [23] . We show here that the coupling with the exter¬ 
nal bath gives rise to damping of the scissors modes and 
in addition that this damping offers the possibility to ob¬ 
serve a dynamical Casimir effect. Moreover, we demon¬ 
strate that the amount of damping also depends on the 
external pumping. 

In this article we study light in nano-fabricated semi¬ 


conductor microcavities. The layout of the paper is as 
follows. In Sec.|ll| we model the semiconductor by one 
conduction and one valence band and describe the in¬ 
teractions between the electrons and holes by an effec¬ 
tive contact interaction, which is especially appropriate 
in the electron-hole plasma regime of interest to us here 
where the Coulomb potential is screened to a short-range 
interaction. We include screening effects by calculating 
the scattering length for the appropriate Yukawa poten¬ 
tial. Hereafter, we consider the coupling between light 
and the semiconductor in Sec. m We show that the cou¬ 
pling results into damping of the light and we demon¬ 
strate that this damping can be characterized by a single 
dimensionless parameter that depends on the external 
pumping. After this, we apply our model and study an 
example of a many-body phenomena of light where the 
advantages of semiconductor microcavities are particu¬ 
larly useful. Namely, in Sec. |IV| we propose to investigate 
the superfluidity of light via the excitation of scissors 
modes. In Sec. lIV A] we consider a sudden rotation of 
the trapping potential of an elongated photon conden¬ 
sate and we calculate the decay rates of the excitations 
in the Thomas-Fermi limit. Hereafter, in Sec. |IVB] we 
propose to measure the density-density correlations of 
the excited light fluid, since the decay products of the 
scissors mode quanta also demonstrate an analog of the 
dynamical Casimir effect. Finally, we end with conclu¬ 
sions and discussion in Sec.jVj 

II. SEMICONDUCTOR 

To investigate Bose-Einstein condensation of light in 
a nano-fabricated semiconductor we in first instance ig¬ 
nore the light and start with a model for a homogeneous 
semiconductor. We consider the following action 


Sscir,^] ( 1 ) 

-^0 J 

(A,a' ^ 

X I4(x- x')(/)ft,,„'(x',r)(^e,a(x,T). 


In this model we only take into account one valence and 
conduction band and i denotes the electron e or hole h 
respectively. The generalization to for instance one con¬ 
duction band and three valance bands is straightforward 
and can be easily achieved once the experimentally rele¬ 
vant semiconductor is known. The electron field and hole 
field are denoted by and (j)* ^ with a and a' repre¬ 
senting the spin degeneracy that is denoted by f or j,. 
The noninteracting Green’s function is defined as 
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where fii is the chemical potential and rrii is the corre¬ 
sponding mass of the electron or hole. Note that the 
band gap of the semiconductor is absorped in our defini¬ 
tions of the chemical potentials. The same is true for the 
bandgap renormalization due to the electron-electron and 
hole-hole Coulomb interactions. The electron-hole inter¬ 
action, however, needs to be considered explicitly due to 
the possibility of exciton formation. 

In a semiconductor there are usually Coulomb inter¬ 
actions. However, since we are mostly interested in the 
highly excited regime where excitons do not exist and 
screening plays a dominant role, we may simplify the 
theory and replace the interaction potential by a contact 
interaction 

I4(x - x')-Vb<5(x - x'), (3) 


A. Interactions 

To obtain a better understanding of the effect of the in¬ 
teractions in our model, we now focus on the many-body 
T-matrix. For simplicity we consider relatively small 
densities of electrons and holes, such that the many- 
body effects, i.e., the Fermi-Dirac distribution function 
in Eq. ([^, can be neglected. In this regime we find 

1 1 rrir^'^ ! ] “ 

^ - ^;7r^ Y- 2 ( 2 : - ep.e+/i + Me + (7) 

where z is the complex center-of-mass energy associated 
with the Matsubara frequencies and nir is the so- 
called reduced mass which is defined as 


with Vo the effective interaction strength that we deter¬ 
mine selfconsistently later on. 

Now we perform a Hubbard-Stratonovich transforma¬ 
tions to the pairing fields and integrate out the fermionic 
fields [21]. Note that we have to perform four differ¬ 
ent transformations and therefore we introduce the fields 
Aq,c('(x,t), of which the averages are given by 


nrir 


memh 


me 


■ ruh 


( 8 ) 


Moreover, ep^e+h is the kinetic energy of the center-of- 
mass motion of the electron and hole involved in the in¬ 
teraction that is given by 


(A„q/(x,t)) = Fo (</>/»,a (x,T)0e.a'(x,T)). (4) 

We apply here the Nozieres-Schmitt-Rink approximation 
and only consider terms up to quadratic order in the pair¬ 
ing fields |25] . which is the simplest approximation that 
correctly incorporates the crossover between an exciton 
gas and an electron-hole plasma that occurs as a func¬ 
tion of excitation, i.e., pumping, of the semiconductor. 
In this approximation we find that the thermodynamic 
potential is given by 

H := Hi + H 2 = ^ ^ In P)) (5) 

P,n 



where are bosonic Matsubara frequencies and Cpy = 
h?p^/2mi is the kinetic energy of the particle or hole. 
The thermodynamic potential thus consists of the sum 
of the ideal electron and hole contributions and a fluctu¬ 
ation correction. The many-body T-matrix in the above 
expression is defined as 

_ I _ ^ J_ — ^ ('C'l 

P)~ Vo Cpv + V./* ^ ’ 

_l_ 1 — jVFD(eP-p',e — Me) — NpD^ep/^h — IJ-h) 

V ihVtji 6p —p'^e T -h /^e M/i 

where Npd{x) = l/(e^^ -I- 1) is the Fermi-Dirac distribu¬ 
tion function. From now on we simplify the notation by 
omitting the arguments of the many-body T-matrix. 


ep,e+h = —( -,-r- (9) 

2(me -I- mh) 

We write the effective interaction strength as Vb = 
mr with a the so-called effective scattering length 
of the potential between the electron and hole. 

From the expression for the many-body T-matrix in 
Eq.0 we obtain two features as a function of the com¬ 
plex energy z. Eirst, we find that for positive values of 
Vo there is a pole on the real axis. This corresponds to 
the exciton that is a bound state of an electron and a 
hole. Second, we find that there is a branch cut on the 
real axis, which corresponds to the continuum of electron- 
hole scattering states. To see this explicitly we calculate 
the density of electrons in the semiconductor. As we 
only consider optical excitations, the electron density Ug 
is equal to the hole density Uh- Hence, we define a carrier 
density n as 

" = ”' = "‘ = -KSr^ 

Since we have an analytical expression for the many-body 
T-matrix, we can also calculate the density analytically. 
By using similar techniques as shown in Ref. |26j , we per¬ 
form the summation over Matsubara frequencies by con¬ 
tour integration. Note that we have to be careful with 
the branch cut and the pole when choosing the contours. 
Ultimately, we find 

n = Uid -f rZex + «sc, (11) 
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FIG. 1: a) Number of electrons in excitons and b) total number of electrons in the excitons and the scattering contributions 
as a function of the inverse of the interaction parameter a for ZnO. The solid and dashed line correspond to a carrier density 
of n = 10^® m“^ and n = 10^'* m“® respectively. We obtain a smooth crossover from Ath/a 0 where the electrons and holes 
behave as an ideal gas and form an electron-hole plasma, to Ath/a ^ 0 where almost all electrons and holes form excitons. 


where 


^ j J dey/eNpuie- (12) 

0(a) /2{me + 


^ex — 9 
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de y/e Nbe{ ^ex); 
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de 
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2/2 -I- (Ath/4aV7r) 


with 0 the Heaviside step function, /3 = l/k-eT^ 

Ath = {2'KfPp/mr) the thermal de Broglie wave¬ 
length, AfBE(a:) = l/(e^ -I- 1) the Bose-Einstein distribu¬ 
tion function and we defined the dimensionless exciton 
energy Cex as 


Eex = e - /3(/ie + yh) 


±_ //wy 

IGtt \ a / 


(13) 


instance and come back to this problem in a moment. 

In Fig. [l] we show the exciton contribution riex and the 
sum with the scattering contribution rise as a function of 
the scattering length. As an example, we take ZnO at 
room temperature and use me = 0.28mo, = 0.59mo 

with Too the bare electron mass, see Ref. [23]. First, we 
note from Eqs. (12) that there are only excitons for posi¬ 
tive values of the scattering length and therefore we only 
show the number of excitons in this regime. For smaller 
values of the scattering length the contribution of the ex¬ 
citons becomes larger. Furthermore, by looking at the 
sum of the exciton and scattering contribution as a func¬ 
tion of the interactions, we note that there is a smooth 
crossover from Ath/a ^ 0 when the electrons behave 
as an ideal gas, to Ath/a ^ 0 where almost all electrons 
form excitons. Therefore, we note that this model for the 
semiconductor describes both the exciton regime and the 
regime where there is a electron-hole plasma. Moreover, 
we note that the sum of the exciton and the scattering 
contribution is an analog for the behavior of the Cooper- 
pair density in the BEC-BCS crossover [27| . 


We see that there are three contributions to the carrier 
density. The first contribution is simply the ideal gas re¬ 
sult. The second part originates from the exciton bound 
state and the third part describes the scattering contin¬ 
uum of electron-hole states. We can calculate the sep¬ 
arate contributions to the density as a function of the 
value of the scattering length. In order to obtain these 
values, we first take a certain carrier density n. Then 
we can find the chemical potentials as a function of the 
scattering length, by using Eqs. 0 and for both 
i = e and i = h. By resubstituting the chemical poten¬ 
tials, we obtain that every contribution to the density 
only depends on the scattering length and the carrier 
density. Therefore, we investigate the effect of changing 
the value of these physical parameters. Of course, be¬ 
cause of screening the effective scattering length will also 
depend on the carrier density, but we ignore this in first 


B. Scattering length 

Up to now we have considered the scattering length as 
a free parameter. However, in a semiconductor the in¬ 
teraction depends on the value of the carrier density due 
to screening effets. Here we only consider the effects of 
static screening that can be described by using a Yukawa 
potential. Therefore, we can make our model selfconsis- 
tent by calculating the scattering length for the Yukawa 
potential. 

The standard procedure for calculating the scattering 
length is via the so-called Born series. However, for a 
Yukawa potential this will lead to divergences and taking 
into account only the first term of this expansion is not 
sufficient. Therefore, we here use a different approach. 
First, we note that the scattering length is associated 
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FIG. 2: The scattering length for ZnO as a function of the 
carrier density by using the Yukawa potential with a screening 
length that is calculated for the ideal-gas parts of the electron 
and hole densities. We obtain that the last resonance occurs 
at the Mott-density riM — 2.3 • 


with a low-energy two-body scattering wavefunction that 
is given by 


^(r) ^ ( 1 - ^ 


(14) 


where A is a normalization constant that is usually taken 
to be 1 for many-body applications. Moreover, if we de¬ 
fine u(r) = riplr) than this function satisfies the following 
radial Schrddinger equation 


h? cPu{r) 
2mr dr'^ 


= V (r)u(r), 


(15) 


where V (r) is a potential and r is the relative coordinate 
between the electron and hole. Note that this is an equa¬ 
tion for the relative motion between the electron and hole 
and therefore we have to use the reduced mass in the 
kinetic energy part. For relatively large distances the po¬ 
tential vanishes and therefore the wavefunction given by 
Eq. (141 is a solution to this equation. Hence, we have 


to solve this equation and look at the solution at large 
distances. To find the scattering length, we use that 


-a = lim 

r—)-oo 


u(r) 

/(r) 


(16) 


Thus by solving the radial Schrddinger equation as writ¬ 
ten in Eq. (15) we can find the scattering length. 


For the semiconductor we assume that there is only 
static screening and therefore we are allowed to take a 
Yukawa potential that is defined as 


assuming that only the unbound carriers screen the po¬ 
tential. Although this approximation leads to an overes¬ 
timation of the screening length, we expect that this is a 
small effect since screening by bound carriers is weaker 
than screening by unbound carriers. 

In this approximation the screening length A, satisfies 
Aj2 = a- 2 + with 


A. 


epCr 

dn 


(18) 


the screening lengths of the electron and hole plasmas 
respectively. Note that in this formula the superscript 
0 of the chemical potential indicates that we calculate 
from the unbound carrier density by treating the elec¬ 
tron and hole as ideal Fermi gases. To find the scattering 
length for the Yukawa potential we still need to specify 
the boundary conditions in Eq. ( |15[ ). Since the potential 
is infinite at r = 0, we must take u(0) = 0. Moreover, we 
set u'(0) = 1. Note that this latter boundary condition 
is simply a matter of normalization and does not affect 
the value of the scattering length. 

In Fig.we display the scattering length as a function 
of the carrier density by determining the screening length 
with the ideal-gas parts of the electron and hole densities. 
For convenience we again consider ZnO with the values 
as stated in Ref. PH ■ We observe that there are multiple 
resonances or equivalently values of the carrier density 
where the scattering length diverges. This is a conse¬ 
quence of the changing value of the screening length for 
different carrier densities. For larger carrier densities, the 
screening length becomes smaller and therefore roughly 
speaking the potential is less deep. Hence, for small car¬ 
rier densities the Yukawa potential supports more bound 
states. Furthermore, every resonance corresponds to a 
bound state and we find that the last resonance occurs at 
the Mott-density um — 2.3-10^^ m“^. After that value of 
the carrier density, the scattering length always remains 
negative. Therefore, this calculation shows that there are 
only excitons for carrier densities that are smaller than 
the Mott-density. Note that this is in agreement with 
previously obtained results and the value for the Mott- 
density is within the range of published data 

For large carrier densities we only have to take into ac¬ 
count the last resonance. Moreover, for carrier densities 
that are close to um we can approximate 


ajn) ^ 2 I 

Ath Um — n 


(19) 


C. Many-body effects 


V{r) 


AireoCr 


c-r/Xs 


(17) 


with e the electron charge, ep the vacuum permittivity, 
Sr the appropriate relative dieletric constant and Ag the 
screening length. We calculate the screening length by 


We are primarly interested in the regime where the 
amount of excitons is negligible, because this is the 
regime where Bose-Einstein condensation of light is pos¬ 
sible. From the results in the last section, we know 
that we therefore have to consider large carrier densities. 
At these high densities the many-body effects, i.e., the 






















6 


Fermi-Dirac distribution functions in Eq. ([^ , become im¬ 
portant. Therefore, we now focus on these contributions 
and we consider the part of the many-body T-matrix that 
has been neglected before, namely 


(27r)3 


dp' 


./VFD(eP-p',e) + NFujip'.h) 

^P —p',e ,h “t“ 


( 20 ) 


In first instance we only consider the imaginary part and 
we neglect the real part. For further purposes we need 
to calculate this many-body T-matrix for —>■ = 

Lo + ie, where e > 0 is infinitesimally small. By switching 
to spherical coordinates and performing the angular inte¬ 
gration, we find that there is only a non-zero imaginary 
part if 


y:= huj + + Hh - ^■p,e+h > 0 . 

The imaginary part is given by 


e(y) 


dpkme iVFD(ep,?i) + (e O h) 


( 21 ) 


( 22 ) 


where P = |P| and (e O h) denotes that there is a similar 
contribution as the first integral where the indices e and 
h are interchanged. Furthermore, 


o'*. 

i^min 


Pn 


= 1 L - 

h y me + mh me + mu 

h V me + mh me + mt 


(23) 


where and p^ax given by similar expression with 
again the indices e and h interchanged. 

By using 


mh ^ 


dpp N-pBiep^h) 


(24) 




we obtain for the full many-body T-matrix, 
•ill 


pMB " ( y) ^^f^i 
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imemh 


X 0(?/)Log 


2 

cosh {huJ + Ce+h}^ + cosh (^^Ce-h 


(25) 


cosh {huJ - Ce+h}'j + cosh (^^Ce-h^ 


where 


„ _ . y/memh - 

Pe+/i — I _ Y y Cp^e-t-^; 


me 


■ mh 


(26) 


Ce—h — pe Ph A~ 


me - mh 


me 


■mh 


[y - ep.e-r/i) ■ 


Note that this results includes the two-body result as 
given by Eq. Q. Moreover, in the context of ultracold 



FIG. 3: The chemical potential of the electrons fie (solid) 
and the chemical potential of the holes fih (dashed) for ZnO 
as a function of the carrier density by using the self-consistent 
scattering length for the Yukawa potential. 


Fermi gases a similar calculation has already been per¬ 
formed in Ref. [5H]. As expected, we agree with that re¬ 
sult if we put me = mh = m and pe = Ph = P- 

Now we use this result to find a self-consistent solu¬ 
tion for the chemical potential of the electrons and holes. 
Recall that 

dV.\ 2 ^ 

= ^y^^FD(ep,0 (27) 

T n 


v{dfLj 




where 




r““(w,p) 


d 


dfi, 


(28) 


As in the previous sections, we only want to consider 
optical excitation. Therefore, we for simplicity from now 
on neglect the dependence of the many-body T-matrix 
on pe — Ph and set 


Ce-h = 


me - mh 


me 


mh 


(y-ep^e+h) ■ 


(29) 


In this case the spectral function is the same for i = e 
and i = h and we again define the carrier density as 

n = rie = n/i. 

The spectral function p(a;,P) = pe(w,P) = ph{oj,P) 
of the pairs contains two different contributions that we 
can separate by considering either ?/ < 0 or j/ > 0. For 
y < 0, we obtain that the Logarithm does not contribute 
in Eq. (251. In this case we can analytically determine 


the spectral function and we find 

y(w,P) 


5 J 

0 (a) . ^th \ 

— S - ep.e+. + Me + y. + j 


( 30 ) 
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By using the result of Eq. (13) we recognize that this 
contribution to the spectral function is simply a delta- 
function at the exciton energy. Again recall that this 
contribution is only present when the scattering length 
is positive and from the results in the previous section 
this implies that for carrier densities larger than the Mott 
density this contribution is not present. The other con¬ 
tribution to the spectral function is obtained by taking 
y > 0. In this case the square root in the expression for 
the many-body T-matrix vanishes and we are left with 
the contribution from the logarithm. This represents the 
scattering continuum of electrons and holes. 

In Fig.j^we show the result for the chemical potentials 
of electrons and holes as a function of the carrier density 
n for ZnO as discussed in Ref. [24j . In comparison with 
the calculation in that article, where the chemical poten¬ 
tials are determined by treating the electron and holes as 
two ideal Fermi gases, we find slightly smaller values for 
both chemical potentials. For example, we find for the 
value of the carrier density that corresponds to popula¬ 
tion inversion, i.e., /ie -I- /r/j > 0, a value that is roughly 
two times larger than stated in Ref. [23]. 


III. LIGHT IN SEMICONDUCTOR 
MICROCAVITIES 


Now that we worked out the model for the semicon¬ 
ductor, we consider the full situation and include the 
coupling to an external light field. Similar to Bose- 
Einstein condensation of light in a dye-filled cavity, we 
can describe the complete dynamics using the Schwinger- 
Keldysh formalism as presented in Ref. |2S] • As will fol¬ 
low from the results in this section, in this case the self¬ 
energy as given in the Langevin field equation in Ref. [29j 
is proportional to the susceptibility of the semiconductor. 
However, also for Bose-Einstein condensation of light in 
nano-fabricated semiconductor microcavities the system 
relaxes towards a steady state that can be described by 
standard equilibrium methods. It is important to realize 
here that throughout our paper we are always considering 
the thermalized case, and in particular are not discussing 
the fully non-equilibrium laser regime of this system. We 
start from the action 


S[ak, Ok, 0k, 0y = 5'sc[0k, 0k 


(31) 


r^P f d 1 

+ d'rak(T) j/i—-I-eph(k) -/Tphj ak(T) 

0Tak(T)0h,p _a(T)0e,k-p,a(T) 


k,p,a ' 


9eh 

vv 


f cl'^«k(T)0*^k-p.aM0)).p,-a('r)- 

1 -'0 
k,p,a 


This action contains several parts. First, 5'sc[0k,0y de¬ 
scribes the semiconductor, which is the momentum-space 
representation of Eq. Q. The second part describes the 



photon energy (eV) 


FIG. 4: The imaginary part of the susceptibility Xo{^) di¬ 
vided by the energy of the photon at the minimum of the 
band Eq as a function of the photon energy for — 
1.7 • 10“®^ J • m^A and Eq = 0.72 eV. The dotted, dashed and 
solid curves corresponds to a carrier density of n = 10^® rn“®, 
n = 5 ■ 10^® m~^ and n = 10^® m“® respectively. 


light with ak('r) and ak(r) the photon fields. Here, eph(k;) 
is the kinetic energy of the photons. By performing a full 
band structure calculation, e.g. see Ref. m, we know 
that 

eph(k)-Fo + ^^fc^, (32) 

with Eq the energy of the photon at the minimum of the 
band with respect to the energy of the band gap of the 
semiconductor, ko the wavenumber at the minimum of 
the band and m the effective mass that arises from the 
local curvature of the band. In equilibrium, the sum of 
the number of electron and photons and the sum of holes 
and photons is conserved. Therefore, in equilibrium we 
have 


Mph — Me 4“ l^h- 


(33) 


The last two terms of Eq. (31) describe the interaction of 


the light with the semiconductor. The first of these terms 
corresponds to the annihilation of electron-hole pairs by 
the electric field and the second term is the creation of 
electron-hole pairs by the electric field. Recall that a 
represents the spin degeneracy of the electrons and holes 
and can be either f or In the interactions terms we 
only consider electrons and holes with opposite spin and 
therefore we only take into account transitions without 
spin-flip. Moreover, we introduced the coupling constant 
geh and we assumed that the coupling is independent of 
energy. Finally, note that the photons and the corre¬ 
sponding momentum k are two dimensional. However, 
the electron and holes are three dimensional and there¬ 
fore p is also a three-dimensional vector. This convention 
will be used throughout the remainder of this paper. 

To obtain the effect of the coupling to the semiconduc¬ 
tor on the behaviour of the photon gas, we use second or¬ 
der perturbation theory in the coupling constant geh- By 












photon energy (eV) photon energy (eV) 


FIG. 5; The spectral function of the photons above and below the Mott density of n ~ 2.3 • 10^^ for g^h — 9 • J • 
and Eo = 0.4 eV. In Figure a) the density is n = 3.1 • 10^^ and Figure b) is for a carrier density of n = 1 • 10^® m“®. Below 
the Mott density there are two different contributions corresponding to the upper and lower exciton-polaritons. However, above 
the Mott density excitons do not exist and therefore there are only pure photons that acquire a finite lifetime through the 
interaction with the electrons and holes. This is the regime in which Bose-Einstein condensation of light is possible. 


following the same steps as in the appendix of Ref. [24] , 
we find the effective action for the photons to be 

S'“®[ak,ak] •= “X! alnhG~^{'k,iuJn)ak,n (34) 

k,n 

= X! “k,n{ - + eph(k) - fj-ph - Xk(w„)|ak,n. 

k,n 


Here the susceptibility Xk(w„) acts as a selfenergy for the 
photons. In the Nozieres-Schmitt-Rink approximation it 
is given by 


Xk(Wn) = 


V 


with 


Xk,p(w„) = Xk.p(Wn) ( 1 - ^ XI XKp'i^n) ) ■ (36) 


V 


p' 


Hence, 

Xk,p(^^n) = 

where 

Xk.p(w„) = 


Xk,p(‘^n) 


1 — jVFD(ek-p,e) — .^FD(ep,^) 

ihuJji Ck—p,e ^p^h T Me 4” 


(35) 


(37) 


(38) 


From the previous section we know the effective in¬ 
teraction strength Vq and the chemical potentials as a 
function of the carrier density. Therefore, we can de¬ 
termine the finite lifetime effects of the photons due to 
the interaction with the electron-hole plasma by calcu¬ 
lating the imaginary part of the susceptibility Xk(^^) for 
every carrier density. In principle we can determine the 


susceptibility for arbitrary photon momentum k. How¬ 
ever, since we are primarily interested in Bose-Einstein 
condensation of the photons, we from now onwards only 
consider k = 0. 

In Fig.j^we display the imaginary part of the suscepti¬ 
bility for geh — 1-7 • 10“^^ J • m^/^. Note that this values 
only changes the absolute value of the imaginary part 
and does not affect the qualitative behaviour. However, 
the actual value for these physical parameter should be 
obtained from experiments on optical spectra of the semi¬ 
conductor. We find that for small carrier densities the 
imaginary part is only nonzero for positive photon ener¬ 
gies. However, for carrier densities that are large enough, 
we also obtain negative contributions for negative ener¬ 
gies. This is a consequence of the fact that only for large 
enough carrier densities there is gain of photons by emis¬ 
sion processes. Moreover, we find that for large carrier 
densities, where there is substantial gain of photons, the 
imaginairy part of the susceptibility is linear for small 
values of the photon energy. Note that this is similar 
to the result for Bose-Einstein condensation of light in 
a dye-filled optical microcavity as presented in Ref. [29] . 
In this regime of large carrier densities we can combine 
the finite lifetime effects of the photons in a single di¬ 
mensionless parameter a, which is simply the slope of 
the imaginairy part of the susceptibility at the origin. 
Note that we considered a homogeneous semiconductor 
coupled to an external light field. However, in the pro¬ 
posed experiment as envisaged here, the semiconductor 
contains airholes, in which there are of course no electron 
or holes. In this case we need to multiply a by 1 — r^air, 
where r^air corresponds to the fraction of the volume of 
the holes of the semiconductor. 

Up to now we only considered carrier densities that are 
larger than the Mott density. In this case we are in the 
pure photon regime, whereas for carrier densities that are 
smaller than the Mott density the photons are coupled to 
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the exciton and therefore exciton-polaritons excitations 
exist. To emphasize the difference between the two den¬ 
sity regimes, we now investigate the spectral function of 
the photons in the two different regimes. We define the 
spectral function of the photons as 


p^(k,a;) = [G'T,(k, a;+)] 


(39) 


where G^(k, ioJn) is defined in Eq. (34) and = uj + ie 
with e infinitesimally small. 

In Fig.we show the dimensionless spectral function 
p-yCkjUj) ■ eV for geh — 9 • 10“^^ J • and Eq = 0.4 eV 
as a function of the photon energy. Below the Mott den¬ 
sity we have two distinct contributions. Due to the cou¬ 
pling of the excitons to the photons, we have two peaks 
corresponding to the lower and upper exciton-polariton 
branches. Note that the photons that are part of the 
upper exciton-polariton acquire a finite lifetime, since 
their energy is larger than the energy threshold for the 
scattering continuum of electrons and holes. However, 
for carrier densities that are larger than the Mott den¬ 
sity, the excitons are no longer present. In this case the 
spectral function only has a pure photon contribution, 
where the interaction with the semiconductor results into 
a finite lifetime of the photon. Therefore, in this spec¬ 
tral function we clearly find the physical differences be¬ 
tween the exciton-polariton and the photon limits of the 
crossover. It is important to note that Bose-Einstein 
condensation of photons appears only in the regime de¬ 
scribed in Fig.j^b. 


IV. SCISSORS MODES 

In the previous section we constructed a model for 
Bose-Einstein condensation of light in nano-fabricated 
semiconductor microcavities. From now onwards we fo¬ 
cus on an example of a many-body phenomenon of light 
that cannot be studied the currents experiments on Bose- 
Einstein condensation of photons, but can be investigated 
in the proposed experimental set-up. Namely, we con¬ 
sider probing the superfluidity of light via the excitation 
of scissor modes. 

We consider a two-dimensional photon gas with effec¬ 
tive mass m and an effective point-like interaction g in 
an external harmonic trapping potential V®’'*(x) with 
trapping frequencies tOx and ujy. Note that in the nano- 
fabricated semiconductor microcavities the harmonic po¬ 
tential arises by systematically increasing the size of the 
holes from the center to the edge of the semiconduc¬ 
tor. Furthermore, in the previous section we showed that 
there is an excitation exchange between the photon gas 
and the electron-hole plasma, but, apart from such pro¬ 
cesses, photons are assumed to be conserved. Therefore, 
we are allowed to consider this quasi-equilibrium photon 
gas in the grand-canonical ensemble and to introduce the 
chemical potential p, of the photons. Moreover, we intro¬ 
duce a dimensionless interaction parameter g that de- 




2;b 


V 










V 


FIG. 6: The cigar-shaped superfiuid photon gas after a sud¬ 
den rotation of the trapping potential by a small angle X- 
The dashed lines indicates the oscillations of the condensate 
and the solid arrows denote the irrotational velocity v of the 
condensed photons. The red wiggly arrows represent the de¬ 
cay products of the scissors mode quanta that can be used to 
demonstrate the dynamical Casimir effect. 


scribes the interactions between the photons. Note that 
an explicit expression for the parameter can be obtained 
from Eq. (31) by applying fourth order perturbation the¬ 
ory in the coupling constant ge,h- However, here we sim¬ 
ply use a phenomenological approach to incorporate the 
interactions between the photons and leave the calcula¬ 
tion of the photon-photon interaction for future work. 

We are interested in the dynamics of a condensate of 
light after a sudden rotation of the trap. In the following, 
we treat the local density and local phase of the conden¬ 
sate separately and we consider a condensate with a large 
number of photons such that we can use the Thomas- 
Fermi approximation. In this approximation the equilib¬ 
rium condensate density no(x) is equal to 


?^o(x) = - 1 — (a;/i?TF,a:) — {v / RtF ,y) 


(40) 


with i “ 2p/muif the Thomas-Fermi radius in the 
corresponding direction. Furthermore, the density is zero 
outside the ellipse that is spanned by these two radii. A 
rotation of the trap by a small angle x results into a 
change in the condensate density 6n(x) that is given by 

5n{x,t) := [no(xO ~ rio(x)] = (7x7/6“*“*, (41) 

where x' are the coordinates after the sudden rotation of 
the trap and C := 2xpiR^p x ~ y)/9- Here we used 

that for an eigenmode the time dependence is harmonic 
with angular frequency w. To obtain more information 
about the fluctuations of the local phase, we consider the 
hydrodynamic equations of the condensate as for example 
can be found in Ref. [30] . By assuming that both the 
velocity of the condensate and the density fluctuations 
(5n(x) are small, we find for the phase of the condensate 


d<(.(x,t) = ^(7xye-“*. (42) 

inu! 

Moreover, for a harmonic trapping potential the fre¬ 
quency of the scissors mode is given by w = . 

In the normal state the situation is different and we 
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have to distinguish between the collisional and collision¬ 
less regime. In the collisional regime the frequency of 
the excitations is the same as in the superfluid phase 
m- However, if the gas is dilute and the interactions 
are weak, the system is in the collisionless regime and 
the analog of the scissors-mode frequency is equal to 
LOx + Wy. By comparing the classical collision rate with 
the trap frequencies, we can make a distinction between 
both regimes. By using the expressions as specified in 
Ref. |19j . we know how to construct the experiment con¬ 
ditions such that we are in the collisionless regime and 
we can distinguish whether the photons are superfluid or 
not by measuring the frequencies of the excitations of the 
photon gas after applying a rotation to the trap. 


A. Damping of scissors modes 


We have just seen that the superfluidity of the pho¬ 
tons can be studied by looking at the frequencies of the 
excitations after applying a sudden rotation to the trap. 
These scissors modes result in oscillations that can di¬ 
rectly be observed in experiments. However, due to the 
coupling of the photons with the electron-hole plasma it 
is worthwhile investigating how this coupling leads to the 
damping of these oscillations. 

We have demonstrated that for a Bose-Einstein con¬ 
densate of light, at large carrier densities the effects of 
the electron-hole plasma can be characterized by a sin¬ 
gle dimensionless damping parameter a that depends on 
the external pumping. Moreover, the photon gas equi¬ 
librates to a steady state that is a dynamical balance 
between particle losses and external pumping. In the fol¬ 
lowing we start from this steady state and we investigate 
the associated decay processes of the scissors modes. In 
particular, we show that the decay rate depend on the 
value of a and thereby on the external pumping. 

To observe properties of the damping, we are primar¬ 
ily interested in configurations that allow for many decay 
processes. For an elongated trap, the difference between 
the energies of two adjacent modes in the direction with 
the small trap frequency is small. Therefore, we expect 
that in that case the energy of the excitations in the 
long direction almost forms a continuum and the scissors 
mode quanta can decay into many other modes. Hence, 
this elongated trap is particularly interesting and we only 
consider this configuration throughout the remainder of 
this paper. A summary of the proposed structure is dis¬ 
played in Fig.[^ 

We are interested in damping processes where fluc¬ 
tuations of the condensate induce the creation of non- 
condensed excitations and therefore we only have to con¬ 
sider the interaction part of the hamiltonian. We substi¬ 
tute for the creation and annihilation operator of the pho¬ 
tons V’(x, t) = ('0(x, t)) -I- (5^/!(x, t) and we only consider 
terms up to quadratic order in the fluctuations 


In this Bogoliubov approximation, we therefore study 

^int(t) = I J (43) 

+ 9 J 


Now we explicitly separate the dynamics of the lo¬ 
cal phase and density of the condensate by writing 
('0(x, t)) = Y^no(x)~+”5n(x(t)e®^'^(’^’*\ with no(x) the 
equilibrium condensate density and both 6n(x, t) and 
t) fluctuations that are known from the calcula¬ 
tions in the previous section. Since these fluctuations are 
small, we only consider the first non-vanishing term in the 
condensate fluctuations. For the harmonic potential that 
is considered here, the parts of the hamiltonian that are 
linear in the condensate fluctuations Sn(x, t) and S4>(x, t) 
vanish, since the fluctuations are odd under y —>■ —y. 
Hence, there is no decay of a single scissors mode quan¬ 
tum at this level of approximation and we have to con¬ 
sider the parts of the hamiltonian that are quadratic in 
the condensate fluctuations. 

We introduce the Bogoliubov amplitudes Uj{x) and 
Vj{x) according to 


t)=Y^ [uj (x)6j (t) - V* (x)S] (t) , 
j 

\^U*(x)b](t) - Vj{x)bj(t) , 

j 


(44) 


where j = 1, 2,... indicates the mode number of the ex¬ 
citation. Note that we neglect the contributions coming 
from j = 0, since we are only interested in the decay into 
non-condensed excitations and the j = 0 term correspond 
to the dynamics of the global phase of the condensate. 
Therefore, we have 


jj'^o 


(45) 




fvjj 




where 


2„,2 


■/ = / dxx y 


3,3 

r2 


u*(x)u*,(x)-u*(x)i;*,(x) , (46) 


Hfj'= / dxx^y^no{x) u*{x)u*,{x) + v^x)v*,(x) 


To make further progress, we need to determine the Bo¬ 
goliubov amplitudes and determine the time dependence 
of the operators bj{t) and bj{t). Hence, we need to solve 
the Bogoliubov-de Gennes equations. 

In two dimensions the general solution of the 
Bogoliubov-de Gennes equations is not known. However, 











11 


in the elongated configuration we expect that the dy¬ 
namics of the long-wavelength photons is approximately 
one dimensional. Therefore, we use the exact solution of 
the Bogoliubov-de Gennes equations in one dimension, 
see e.g. Ref. m, in order to make a proper variational 
ansatz for the form of the non-condensed fluctuations 
in the elongated conhguration. Thus, we consider the 
ansatz 


Uj{2.) = 


2^ 

Vj(z) = -^ 

' 2^ ^ 




yr^- 


13. 


- |z 

/3. 




' Pjix), (47) 

Pj{x), 


with z = {x,y) = {x/RTF,x,y/RTF,y) and Pj{x) the j- 
th Legendre polynomial. Furthermore, the constants Uj 
and 13j are given by 

3 yj ~RTF,xRTF,y 

o. _ _ \ _ j ^^3 

^ 2y^RTF,xRTF,y Y 

with Qj = u}x\/j(.j + l)/2 the energy eigenvalues of the 
one dimensional problem. Also Rtf, a; and RxF.y corre¬ 
spond to the Thomas-Fermi radius of the condensate in 
the specified direction. Moreover, we defined the con¬ 
stant Cj which is given by 

Cj = J ^dx{l-x^y/'^ [Pj{x)]^ . (49) 


This ansatz corresponds to a solution for an elongated 
trap where the excitation only propagates in the direc¬ 
tion with the small trapping frequency. Note that if from 
the start we would have removed the j/-dependence, our 
ansatz simplifies to the exact solution in one dimension 
as given in Ref. m- 

Since the photons are in a good first approxima¬ 
tion equivalent to non-relativistic particles with an ef¬ 
fective mass m and point-like interaction with strength 
g, we consider the following action within the functional- 
integral formalism in the Bogoliubov approximation 


S[(l}*,(j3]= J J dx(j)*{: s.,t)G ^(x, r)(()(x, r) (50) 
+ I / J dxno(x) |[(/)*(x,r)]2-h [(/)(x,t)]2| . 

Here we defined 

G"^(x,r) = h— - -f H®’^*(x) 2gn(x). (51) 

Furthermore, (/>(x, r) and 0*(x,t) are the fields that de¬ 
scribe the fluctuations that originate from the Bogliubov 
approximation ^/>(x,r) = ('!/'(x, r)) -|-(^(x,t), no(x) is the 



j 


FIG. 7: The decay rate of two scissors mode quanta into a 
non-condensed mode Wj = + Sjji) with mode 

number j for y — 10 LOy, cuy = WlOx and a = 10“^. 


condensate density, u(x) equals the total density and 
yext(x) corresponds to the external trapping potential. 
We consider the usual Bogoliubov transformation with 


<(>(x,f) = \uj{x)by{T) - u*(x)6*(t) 


(52) 


and the ansatz from Eq. (47) for the coherence factors. 
This allows us to rewrite 




+ h 


dr 




Gjj',n Gj 

1 J 


_ Gi,i'.i2 Gj 


0 

). 

’ bj'A) ’ 

0 ~Gjj' _ 


[W\ 


(53) 


with 


G,,,,=/<ixh(xK,(x)-„,(xK,(x)]. (.54) 

Gjj'.ii = J c^xno(x) [uj(x)u^/(x) -buj/(x)z;j(x)] 

+ /dx[.,(x)G-(x)„,,(x) + .,(x)G-(x)»,(x) 

= g y dxno(x) [uj{x)uf{x) -|-Uj(x)uj-(x)] 

— / dx Uj (x)G“^(x)uj/(x)-|-Uj(x)G“^(x)uj'(x) 


Here we used that Un{'^) and t'n(x) are real. Further¬ 
more, 


G ^(x)=G ^{x,T)-n^ = -^-^+gno{x), (55) 

where we used the Gross-Pitaevski equation in the 
Thomas-Fermi limit. Moreover, in this limit we can also 
neglect the second-order derivative with respect to y. 
This is valid as G~^ acts on Un(x) and v„(x) that only 
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have y-dependence in the density. 

We consider Fj{x) = Uj{x) + Vj(x) and Gj(x) = 
Uj(x) — Vj(x) with Uy(x) and Vj(x) as given by Eq. (47). 
Then, 






4m 


1 — a;^ 


G,(x), 


where we used the results of Ref. [26] and the fact that 
we are in the Thomas-Fermi limit and therefore neglected 
the derivatives of the densities. Namely, all these deriva¬ 
tives result into second-order derivatives of densities or 
terms that contain nQ(x)/no(x), which are all small in 
the Thomas-Fermi limit. Similarly, we find that 


^2 92 

2m dx^ 




1 - x2 - y2 Wj 

l-£2 2 


Now, we define 

J ^dxPj{x)Pj>{x)Vl-x^. 


Then, 


(57) 


(58) 


for every j = 0,1,.... Hence, 

bj (r) = c;e- + dj , (62) 

6*(r) = 

Thus the energy of the excitations for j = 0,1,2,... is 
equal to 


hujj = hflj a/2/3 = hujx a/jXj -I- 1)/3. (63) 


To quantify the effect of neglecting the non-diagonal 
terms, we numerically determine the eigenvalues of the 
full problem, where we consider terms from y = 1 to 
J = 75. The largest error of roughly 15 percent occurs 
for J = 1. For the other values of j the error is even 
smaller and less than 10 percent. Therefore, it is indeed 
a good first approximation to neglect the non-diagonal 
terms. Furthermore, by resubstituting the expressions in 
Eqs. (62) into Eqs. (611 we obtain 


d* = - ( 5 -f 2v^) c*, (64) 

Cj = ^2-\/6 — 5^ dj. 


G, 


hi' 


1 

2 


fif + 1 ) , 4 / JU + 1 ) 

j{j +1) V 


(59) 


Gj.i ',11 


Mlj! 

~ir 



fif +1 ), 


hrij 


Mlji 

^ir 


fif + 1), 
j(j + i) 


For now we ignore the coupling between modes where 
j ^ f. Note that this is in general a good approximation 
as the integral Ijji has a maximal value for j = j' and 
gradually decreases if j is further and further away from 
j'. This allows us to rewrite 



1 

2 


j ^ 

drY^bjir) (h-^ 
j ^ 


5h^ij \ 

) 

bMlj \ 

J 


+ ^/ dT^{bj{T)f P {b*{T)f 
3 


bj{T) 


(60) 


To complete our variational approach we determine the 
equations of motion, which allows us to determine the 
time dependence of the operators bj^r) and 5*(t). We 
find for the equations of motion 


9 


5hQj \ 

^ ) 

bhVtj \ 

J 


bAr) + ^b){r) = Q, 

bj{r) - ^bj{T) = 0, 


(61) 


Hence, up to an overal normalization constant 


where 


(65) 


(66) 


Uj(x) = Uj(x) + (5 — 2-\/6^ v*(x), 
v*{x) = v*{x) + (5 — 2^6^ Uj{x). 


By using these expressions for the Bogoliubov ampli¬ 
tudes, we find for Hjji as defined in Eqs. (45) and (46), 


H - 

mtu 


2,,2 ( 7 ?XF,y '^TF,; 



p2 p2 

^TF,y^TF,x 

wv 2uj 


(67) 


ry2,3 


7O.I yO,l 


16/^2 


72.7 

'‘ 3 , 3 ' 


35 fko^yhujjfujjj I 

where we used the shorthand notation 

/ I / \’t./2 

da;i™^l — Pj{x)Pj>{x). (68) 


ryin^ri _ 


To investigate the damping of the scissors modes, we 
are interested in the transition rate for creating two ex¬ 
citations with frequency LOj and ujji from the vacuum 
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through the decay of two scissors mode quanta. By ap¬ 
plying Fermi’s Golden Rule, we obtain 


use Eq. (44) to rewrite the result in terms of the Bogoli- 
ubov amplitudes w(x) and u(x). We find 




^(1 + 


p{ujj +UJf), 


(69) 


where 


p{u}j+U}f) 


i_ +ajy) _ 

ttH {ujj + ujji — 2w)^ -I- [a(ujj + 


(70) 


Here we introduced a final density of states p{E) to in¬ 
corporate that, due to the interaction with the molecules, 
there is a probability that the photon is in a state with 
an energy that is within a small band around LOj + LOj/. 
Note that for the current experiments we have that 
-I- hujji) 1 with P = l/feelE the inverse of the 
thermal energy. Therefore, we used the low-energy ap¬ 
proximation for the density of states [5^ . As the dimen¬ 
sionless damping parameter a is small, we directly satisfy 
Wjj' ~ 0 to a very good approximation if the photons 
scatter into a state with energy outside the small band 
around ujj +u}ji. Note that the nonzero value of a makes 
our calculation specific to dissipative Bose-Einstein con¬ 
densates and not immediately applicable to a cold atomic 
gas, where a = 0 and Beliaev damping of the scissors 
modes is only possible in the presence of fine-tuned de¬ 
generacies [ 21 ]. 

In Fig.we show the decay rate Wj = J2j' + 

Sjj') for p = 10 ujy, LUy = 10 ujx and a = 10“^. Because 
the frequency of the scissors mode w is roughly equal to 
Wi 7 , we obtain a peak for j = 17. Furthermore, we find 
that the decay of the scissors modes indeed leads to the 
population of several non-condensed modes. 


B. Density-density correlation function 

We now consider the situation that the scissors modes 
are being excited and we consider the density-density cor¬ 
relation function in the operator formalism. Thus, we 
define 


5(2)(x,x',t) = 


(p(x,t)p(x',t)) 

(/5(x,t))(/3(x',t)) 


- 1 , 


(71) 


where p(x, t) = is the density operator. 

We again take ■i/'(x,<) = {'ip{x,t)) + and we ex¬ 

plicitly separate the fluctuations that are described as 
scissors modes and Bogoliubov excitations by writing 
i5i/)(x, t) = i5i/'s(x,t) -I- 5 '(/;b(x, t). As a consequence, the 
density-density correlation function contains the density- 
density correlations from the scissors modes and also a 
term from the density-density correlations between the 
Bogoliubov modes. From now onwards we take y = y' = 
0 such that the contribution of the scissors modes van¬ 


ishes, as can be seen explicitly in Eqs. (41) and (42). 


For the part with the Bogoliubov excitations we can 


g(2)(x,x',t) = -^====^(m„(x)-u„(x)) (72) 
Vno(x)no(x') ^ 

X (m„(x') - r;„(x')) 1 1 -f 2 /S]((t) 6 „(t)^ 


where Nn{t) = is the number of excitations 

in a mode n at time t. The contribution of 5 *-^^ (x, x', t) 
that is independent of Nn{t) has an ultraviolet divergence 
that can be resolved by an appropriate subtraction [32] . 
However, generally this part is neglible compared to the 
contribution that depends on Nn{t) as the photons are 
at room temperature. Thus, in the following we neglect 
this so-called quantum contribution. 

To find the number of excitations Nj{t), we solve the 
following coupled system of equations 


dN,{t) 

dt 



dN,{t) 

dt 


= 0 , 


(73) 


dN,{t) 

dt 


3' 


where N^it) denotes the number of scissors mode quanta 
and Wjj! is the decay rate that is defined in Eq. (69). As 
a lower limit we neglect the Bose stimulation factors in 
the rate equations and in this approximation we find 


m) 




(74) 


Now we use the decay rates as displayed in Fig.and as 
an illustration we consider a time t such that Nj{t) = 
WjlWjnax- The corresponding contribution of the decay 
of the scissors modes to the density-density correlation 
function is displayed in Fig.[^ Note that in experiments 
there is always another contribution from the thermal 
background, which is given by replacing (bl^{t)bn{i)) by 
the Bose-Einstein distribution function at energy 
Therefore, to obtain the result of this figure experi¬ 
mentally, the contribution from the thermal background 
should be subtracted. 

This is only possible when the contribution originat¬ 
ing from this decay process is large enough compared 
to the background contribution. If we consider a time t 
such that Nj{t) = for x x' the value of the 

background contribution is at least two orders of magni¬ 
tude larger and therefore we expect that in this region 
it is difficult to observe. However, for \x — x'\ > 0.2 the 
largest value of the thermal background is maximally a 
factor of ten larger than the scissors mode contribution. 
The exact time scale at which the number of excitations 
Ni (t) takes this value depends on many parameters such 
as the rotation angle y, the number of scissors modes 
and the trapping frequencies. Since in the current exper¬ 
iments the trapping frequencies are of the order of GHz, 
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we obtain that this condition can be fulfilled within the 
nanosecond regime. 

We calculated the number of excitations Nn{t) while 
ignoring the Bose stimulation factors in the rate equa¬ 
tions. To compare with experimental results, the incor¬ 
poration of these additional terms can be important. In 
first approximation these factors can be incorporated by 
replacing the number of excitations by their expectation 
value. As this renormalization increases the decay rate, 
we obtain that within the nanosecond regime the signal 
of the decay of the scissors modes is comparable to the 
background for \x — x'\ > 0.2. Hence, we expect that 
the effect of the decay of the scissors modes can be dis¬ 
tinguished from the background and is observable in the 
density-density correlation function. 

In the elongated configuration, as envisaged here, 
the Bogoliubov amplitudes are proportional to Legen¬ 
dre polynomials. Since these are approximately stand¬ 
ing waves, the scissors mode decays predominently into 
a pair of excitations consisting of an excitation with a cer¬ 
tain local momentum k and —k. This then also explains 
the correspondence with the dynamical Casimir effect as 
an external perturbation, in this case a sudden rotation 
of the trap, leads to the creation of phonon pairs from 
the vacuum. Therefore, a measurement of this density- 
density correlation function would give a demonstration 
of an analog of the dynamical Casimir effect in a Bose- 
Einstein condensate of light, which up to now only is 
considered in atomic and exciton-polariton condensates 

[33H35]. 


V. DISCUSSION AND CONCLUSION 

In this work we discussed a model for a semiconductor 
that qualitatively contains the correct crossover physics. 
However, to obtain the correct quantitative result, several 
improvements have to be made. First, in our model we 
only take into account one conduction and valence band. 
In a realistic semiconductor there are multiple bands that 
all have to be included. Moreover, we have not taken 
into account the band-gap renormalization when self- 
consistently determining the chemical potential. How¬ 
ever, for quantitative agreement it will be important to be 
more careful about this and include the effect of Coulomb 
screening on the renormalization of the band gap. 

Another simplification is the use of the contact interac¬ 
tion for the interactions between the electrons and holes. 
Normally, these interactions are described by a Yukawa 
potential due to screening effects in the semiconductor. 
Although we use the scattering length for the Yukawa 
potential as input for the strength of the contact inter¬ 
action, taking a contact interaction for the interactions 
between the electrons and holes is still a rough approxi¬ 
mation at low carrier densities where screening is not very 
effective. However, we are primarily interested in large 
carrier densities and in this regime the Yukawa potential 
behaves more and more like a contact interaction. We 



X 


FIG. 8: The density-density correlation function of the Bo¬ 
goliubov excitations 8no^^ x', t), where I2mg 

with ^ the so-called coherence length and no the condensate 
density in the center of the trap, in terms of {Yux^'f'hwy/ 
for y = y' = 0. 


verified that in this regime the susceptibility for Yukawa 
interactions and the contact interaction agree quite well. 
For a quantitative agreement at all carrier densities, the 
Yukawa potential has to be taken into account. Another 
extension of our model is to consider dynamical screen¬ 
ing effects, which become important at very high carrier 
densities [5B] . 

In conclusion, we considered Bose-Einstein condensa¬ 
tion of light in nano-fabricated semiconductor microcav¬ 
ities. We modeled the semiconductor as a two-band sys¬ 
tem consisting of electrons and holes that interact via 
a contact interaction. To incorporate screening effects, 
we use the scattering length for the Yukawa potential 
as input parameter for the strength of the contact in¬ 
teraction. We demonstrated that this model contains a 
qualitative description of the regime with and without ex- 
citons. Moreover, we have shown that if we couple light 
to the semiconductor, for large carrier densities the finite 
lifetime effects of the photons can be characterized by a 
single dimensionless parameter a, which is proportional 
to the slope of the imaginary part of the susceptibility 
at zero energy. Hereafter, we have proposed to probe 
the superfluidity of the light in the nano-fabricated semi¬ 
conductor microcavities via the excitation that the scis¬ 
sors modes. By using Fermi’s Golden Rule and a vari¬ 
ational ansatz to calculate the Bogoliubov amplitudes, 
we determined the decay rates of the scissors modes into 
the non-condensed excitations. Finally, we have demon¬ 
strated that the density-density correlation function of 
the excited light fluid exhibits an analog of the dynami¬ 
cal Casimir effect. 
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